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Introduction 


The  investigation  of  heat  and  mass  transfer  proeess  around  a  single  hot  partiele 
surrounded  by  vapor  layer  and  liquid  is  aetual  due  to  a  study  of  vapor  explosion  in  the  time  of 
hard  aeeidents  on  the  nuelear  reaetors  and  metallurgy  industry.  The  proeess  of  heat  and  mass 
transfer  of  vapor  bubble  at  surrounding  liquid  was  investigated  well  [1].  However,  insuffieient 
attention  was  devoted  a  similarly  proeess  in  the  ease  of  presenee  of  a  hot  partiele  in  the  bubble. 
The  models  deseribing  this  proeess  were  proposed  in  the  works  [2,  3  and  4].  Unfortunately  the 
questionable  assumptions  were  used  in  the  papers  [2,  3]  that  ean  restrieted  the  field  of 
applieations  proposed  ones.  The  model  from  the  work  [4]  was  developed  just  for  the  some 
eases. 

The  aim  of  this  work  is  to  fill  up  the  mentioned  gaps. 


The  basic  equations 

Consider  the  system  of  equations  of  model  that  deseribe  the  heat  and  mass  transfer  around 
hot  partiele  surrounded  by  vapor  layer  and  liquid.  The  equations  of  heat  influx  to  a  partiele, 
vapor  and  liquid  ean  be  written  as 
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where  p° ,  7)  (i  =  d,g,l)  are  the  densities  and  temperatures  of  partiele,  vapor  and  liquid;  Wj  , 
cj,  Xj  {j  =  d,l)  are  the  radial  veloeities,  the  speeifie  heats  and  the  heat  eonduetion 
eoeffieients  of  partiele  and  liquid;  Cpg  is  the  adiabatie  speeifie  heat  of  vapor;  d  and  R{t)  are 

the  radii  of  partiele  and  vapor  bubble  aeeordingly.  The  liquid  and  partiele  eonsider 
ineompressible.  It  is  supposed  that  the  pressure  in  vapor  pg  is  homogeneous  in  spaee  and  the 
dependenee  of  thermal  eonduetivity  on  temperature  is  linear: 


The  vapor  mass  balanee  equation  is 
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In  the  case  of  incompressible  liquid  we  have: 
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where  is  the  radial  velocity  of  liquid  near  to  a  bubble  surface. 
We  use  the  ideal  gas  state  equation  for  vapor  in  a  bubble: 
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The  Clapeyron  -  Clausius  equation  relates  the  pressure  in  a  vapor  and  the  saturation 
temperature  : 
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In  this  equations  y  ,  I  and  p°g^  are  vapor  adiabatic  exponent,  specific  heat  of  vaporization  and 

vapor  density  on  the  saturation  line  accordingly. 

The  Rayleigh  -  Lamb  equation  describes  the  oscillations  of  bubble: 
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where  V/  and  a  are  kinematic  viscosity  of  liquid  and  surface  tension  coefficient,  is 


pressure  in  liquid  far  from  bubble,  j  and  are  intensity  of  phase  transition  and  velocity  of 
vapor  on  boundary  surface. 

The  system  of  equations  (1)  -  (9)  is  complete.  The  initial  and  boundary  conditions  are 
required  to  set  for  posing  a  problem.  The  boundary  conditions  for  the  equations  (1)  -  (3)  are: 


r  -  0  : 

dr 

(10) 

r  =  d  : 

^  dr  dr  ’ 

II 

(11) 

r  =  Rit): 

;  57)  dTg  _ 

Pg  =Pl  =PsiPg)’ 

(12) 

r  =  CO  : 

dr 

(13) 

The  boundary  conditions  (12)  on  bubble  surface  describe  the  phase  transition  liquid-vapor  in  a 
quasi -equilibrium  approximation.  Therefore  first  boundary  condition  in  (12)  is  considered  as  the 
equation  for  calculation  of  intensity  of  phase  change  j  .  There  is  the  following  requirement  for 

velocity  of  vapor  Wg  : 

r  -  d  \  Wg  =  0 .  (14) 

It  is  supported  that  initial  temperature  in  a  particle,  vapor  and  liquid  is  linear  function. 
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Transformation  of  eqnations  system 

Similar  to  [1]  we  transform  the  system  of  equation  (1)  -  (14)  to  a  convenient  form  for 
numerical  simulations.  First  of  all  we  note  that  the  equation  (6)  is  inconvenient  for  a  calculation 

of  pressure  in  a  vapor  because  the  density  yc>°  and  the  temperature  Tg  is  depended  on  spatial 
coordinate  r  .  Besides  we  have  not  explicit  equation  for  the  calculation  of  vapor  velocity  Wg  . 
The  equations  for  Wg  and  Pg  can  be  obtained  by  the  substitute  of  expressed  from  equation  (6) 
vapor  density  in  the  equation  (4): 
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Then  the  expression  (15)  is  substituted  in  the  equation  of  heat  influx  to  the  vapor  (2).  The 
retrieved  expression  is  integrated  on  r  from  d  up  to  R{t) .  The  result  of  these  transformations 
is: 
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Repeat  again  the  integration  in  limits  from  d  and  up  to  r'  <  R{t)  and  account  the  condition 
(14)  we  obtain  the  expression  for  the  radial  velocity  of  vapor: 
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Flere  the  prime  over  space  coordinate  r  is  neglected. 

We  use  new  system  of  coordinates  where  the  surface  of  bubble  is  fixed  by  following 
transforming  operators: 
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The  system  of  equations  (1)  -  (9)  and  equations  (16)  -  (17)  for  Wg  and  Pg  may  be 
written  in  the  new  coordinates  as 
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In  these  equations  the  time  z  is  renamed  as  t  and  the  temperature  conduction  coefficients  of 
particle,  vapor  and  liquid  is 
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The  boundary  conditions  (10)  -  (14)  is  transformed  as 
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7  =  1  (1  =  0): 
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We  use  seeond  order  approximation  of  derivatives  on  spaee  for  the  system  of  equations 
(20)  -  (28)  with  boundary  eonditions  (29)  -  (32).  For  integration  on  time  we  use  the  explieit 
differenee  seheme  based  on  Runge  -  Kutta  algorithm  with  the  forth  order  of  aeeuraey.  This 
seheme  was  proposed  in  work  [5]  at  the  solution  of  problem  of  heat  and  mass  transfer  of  vapor 
bubble  in  liquid. 


Conclusion 


The  model  is  proposed  that  deseribes  the  heat  and  mass  transfer  around  hot  partiele 
surrounded  by  a  vapor  layer  and  liquid.  It  was  supposed  that  the  dependenee  of  thermal 
eonduetivity  in  a  vapor  on  temperature  is  linear,  the  pressure  in  vapor  is  homogeneous  in  spaee. 
The  liquid-vapor  phase  transition  was  taken  into  aeeount  in  a  quasi-equilibrium  approximation. 
The  liquid  and  the  hot  partiele  were  eonsidered  ineompressible. 

The  proeedure  of  numerieal  integration  of  system  of  equations  is  developed  that  is  based 
on  the  explieit  differenee  seheme  with  the  forth  order  of  aeeuraey  on  time  and  seeond  order 
aeeuraey  on  spaee. 
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